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Abstract 

We reevaluate the hadronic contribution to the running of the QED fine structure constant a(s) 
at s = M|. We use data from e + e~ annihilation and r decays at low energy and at the qq thresh- 
olds, where resonances occur. Using so-called spectral moments and the Operator Product Ex- 
pansion (OPE), it is shown that a reliable theoretical prediction of the hadronic production rate 
R(s) is available at relatively low energies. Its application improves significantly the precision on 
the hadronic vacuum polarization contribution. We obtain Aahad(Af|) = (277.8 ± 2.6) x 10~ 4 
yielding a~ 1 {Ml) = 128.923 ± 0.036. Inserting this value in a global electroweak fit using 
current experimental input, we constrain the mass of the Standard Model Higgs boson to be 
M ffig gs = 129±<^ 3 GeV. Analogously, we improve the precision of the hadronic contribution to 
the anomalous magnetic moment of the muon for which we obtain a£ ad = (695.1 ± 7.5) x 10~ 10 . 
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1 Introduction 



The running of the QED fine structure constant a(s) and the anomalous magnetic moment 
of the muon are famous observables whose theoretical precisions are limited by second 
order loop effects from hadronic vacuum polarization. Both magnitudes are related via 
dispersion relations to the hadronic production rate in e + e~ annihilation, 

of \ ^tot(e + e--> hadrons) 3s + _ 

R(s) = — — r— = yCTtot{e^e -> hadrons) . (1) 

cr (e + e _ — > fi + ^~) Ana* 

While far from quark thresholds and at sufficiently high energy R{s) can be predicted 
by perturbative QCD, theory may fail when resonances occur, i.e., local quark-hadron 
duality is broken. Fortunately, one can circumvent this drawback by using e + e~ annihi- 
lation data of R(s) and, as recently proposed in Ref. hadronic r decays benefitting 
from the largely conserved vector current (CVC). 

There is a strong interest in the electroweak phenomenology to reduce the uncertainty 
in a(M|) which at present is a serious limit to further progress in the determination of 
the Higgs mass from radiative corrections in the Standard Model. The most constraining 
observable so far has been sin 2 0y^ obtained from leptonic asymmetries at the Z pole with 
an achieved precision of (Asin 2 Ojy) C xp = 0.00022 ||. The uncertainty from the currently 
used a(M|) value translates into (Asm 2 Q w ) a = 0.00023 0, justifying the present work. 

In this paper, we extend the use of the theoretical QCD prediction of R(s) to energies 
down to 1.8 GeV. The reliability of this approach is justified by applying the Wilson 
Operator Product Expansion (OPE) |4j] (also called SVZ approach ||) and fitting the 
dominant nonperturbative power terms directly to the data by means of spectral mo- 
ments. An analogous approach has been successfully applied to the theoretical prediction 
of the t hadronic width, R T , at M T ~ 1.8 GeV in order to measure the strong coupling 
constant a s (M 2 ) [§ ft |, |. 

The paper is organized as follows: first a brief overview over the formulae used is 
given, then the spectral moments are defined and evaluated experimentally, and the cor- 
responding theoretical evaluation is fitted. On this basis, theoretical predictions of the 
hadronic vacuum polarization contribution to the running of a(M| ) and to the anomalous 
magnetic moment of the muon, (g — 2)^, from various energy regimes are determined and, 
with the addition of experimental data, final values for a(M|) and = (g — 2)^/2 are 
determined. 



2 Running of the QED Fine Structure Constant 

The running of the electromagnetic fine structure constant a(s) is governed by the renor- 
malized vacuum polarization function, Il 7 (s). For the spin 1 photon, n 7 (s) is given 
by the Fourier transform of the time-ordered product of the electromagnetic currents 
j£ m (s) in the vacuum {q»q v - q 2 g^)Ii 1 (q 2 ) = i J d A x e iqx (0\T(j» m (x)j^ m (0))\0) . With 
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Aa(s) — 47raRe [II 7 (s) — II 7 (0)], one has 

a(s) = 



a(Q) 



Aa(s) 



(2) 



where 47ra(0) is the square of the electron charge in the long- wavelength Thomson limit. 
The contribution Aa(s) can naturally be subdivided in a leptonic and a hadronic part. 



The leading order leptonic contribution is given by 



A aiop (Ml) = ^E 
on „ 



me 



Mi, 



314.2 x 10' 



(3) 



Using analyticity and unitarity, the dispersion integral for the contribution from the light 
quark hadronic vacuum polarization Aahad(-^z) reads 



Ml 



Air 2 a 



Re / ds 



0"had(S 



Ami 



Ml 



it 



(4) 



where a^^s) = IGir a (s)/s ■ ImII 7 (s) is given by the optical theorem, and ImIT 7 stands 
for the absorptive part of the hadronic vacuum polarization correlator. Through Eq. (JED, 
the above dispersion relation can be expressed as a function of R(s): 



Aa had (Mi) 



iraMl 



Re / ds 



R(s) 



Ami 



s(s - Ml) - ie 



(5) 



Employing the identity l/(x' — x — ie) e ^ = P{l/(x > — x)} + in5(x' — x), the integrals (f|) 
and fl5|) are evaluated using the principle value integration technique. 



3 Muon Magnetic Anomaly 

It is convenient to separate the prediction a^ M from the Standard Model into its different 
contributions 

af = af D + ajr d + , (6) 

where aJ ED = (11 658 470.6 ±0.2) xl0~ 10 is the pure electromagnetic contribution (see [O] 
and references therein), a^ ad is the contribution from hadronic vacuum polarization, and 
a^ eak = (15.1 ± 0.4) x 10~ 10 |11], [L2|, [T3| accounts for corrections due to the exchange of 



the weak interacting bosons up to two loops. 

Equivalently to Aa^diMl), by virtue of the analyticity of the vacuum polarization 
correlator, the contribution of the hadronic vacuum polarization to can be calculated 



via the dispersion integral [14 



ttT = m i ds mm. (7) 



Ami 
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Here K(s) denotes the QED kernel 



K(s) = x 2 (l - y J + (1 + x) 2 (l + (m(l + x) - x + y j + j^y^ 2 lnx (8) 

with a; = (1 — + and /3 = (1 — 4m 2 /s) 1//2 . The function if(s) decreases 

monotonically with increasing s. It gives a strong weight to the low energy part of the 
integral (f|). About 91% of the total contribution to a^ ad is accumulated at cm. energies 
y^s below 2.1 GeV while 72% of ajj ad is covered by the two-pion final state which is 
dominated by the p(770) resonance. Data from vector hadronic r decays published by the 
ALEPH Collaboration provide a very precise spectrum of the two-pion final state as well 
as new input for the lesser known four-pion final states. This new information improves 
significantly the precision of the a^ ad determination [p]]. 



Theoretical Prediction of R(s 



The optical theorem relates the total hadronic width at a given energy-squared Sq to the 
absorptive part of the photon vacuum polarization correlator 

R{s ) = 12vrlmn(s + ie) . (9) 

Perturbative QCD predictions up to next-to-next-to leading order are available for the 
Adler D-function JL6J which is the logarithmic derivative of the correlator II, carrying all 
physical information: 

D(s) = -12tt 2 S ^P . (10) 
ds 



This yields the relation 



\s\=so 



where the contour integral runs counter-clockwise around the circle from s = Sq — ie to 
s = Sq + ie. Choosing the renormalization scale to be the physical scale s, additional 
logarithms in the perturbative expansion of D are absorbed into the running coupling 
constant a s (s). The (massless) NNLO perturbative prediction of D reads then |L7] 
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D P (-s) = N c Y / Qf 
f 




(12) 



where Nc = 3 for SU (3)c and Qf is the charge of the quark /. The coefficients are d = 1, 
di = 1.9857-0.1153 n f , d 2 = d 2 +Pi7r 2 /4:8, (3 = ll-2n//3 and d 2 = -6.6368-1.2001 n f - 
0.0052 n 2 f - 1.2395 (E/ Qf) 2 /N c Ef Q) with n f the number of involved quark flavours^. 
The running of the strong coupling constant a s (s) is governed by the renormalization 

3 The negative energy-squared in Dp(—s) of Eq. ( fL2| ) is introduced when continuing the Adler function 
from the spacelike Euclidean space, where it was originally defined, to the timelike Minkowski space by 
virtue of its analyticity property. 
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group equation (RGE), known precisely to four-loop level ]TS| . 



Using the above formalism, R(sq) is easily obtained by evaluating numerically the 
contour-integral ([y]). The solution is called contour- improved fixed-order perturbation 
theory (FOPT CI ) in the following. Another approach, usually chosen, is to expand a s (s) 
in Eq. (|12|) in powers of a s (s ) with coefficients that are polynomials in ln(s/s ): 



CMS 



71 



7 A) hi— 

7T 4 s 




^iln- 

8 s 



' a s (s) 



71 



4^ In- - In 2 - + ^ 3 In 3 - 

128 s 64 s 64 s 




(13) 



Inserting the above series with the Dp-function in Eq (|TT|) and keeping terms up to oder 
leads to the expression 



R(s ) 



f 



71 



71 



+ d 2 



a s {s 



71 




O (afro)) 



(14) 



where the difference to Dp is of order only. The solution (Tjj ) will be referred to as 
fixed-order perturbation theory (FOPT). There is an intrinsic ambiguity between FOPT 



and FOPTci- The numerical solution of the contour-integral (11) involves the complete 
(known) RGE and provides thus a resummation of all known higher order logarithmic 
terms of the expansion ( |I~3"D (see Ref. |19j for comparison). Unfortunately, it is unclear if 
the resummation does not give rise to a bias of the final result. 



Quark Mass Corrections 

Quark mass corrections in leading order are suppressed as ~ m^(s)/s, i.e., they are small 
sufficiently far away from threshold. Complete formulae for the perturbative prediction 




containing quark masses are provided in Refs. |2lJ [H], 22] for the correlator II(s) and R(s) 
up to order a s exactly and numerically using Pade approximants to order a 2 s . We will 
learn from the numerical analysis that it suffices for the required level of accuracy to use 
the following expansion as an additive correction to the Adler D-function [23] 

iW-«) = -^£Qj^( 6 + 28 ^ + (294.8 - 12.3^)1^] ] . (15) 

The running of the quark mass rrif(s) is obtained from the renormalization group and is 
known to four- loop level f2"4]| . 

When using perturbative QCD at low energy scales one has to worry whether contri- 
butions from nonperturbative QCD could give rise to large corrections. The break down 
of asymptotic freedom is signalled by the emergence of power corrections due to nonper- 
turbative effects in the QCD vacuum. These are introduced via non-vanishing vacuum 
expectation values originating from quark and gluon condensation. It is convenient to use 
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the Operator Product Expansion (OPE) |5], [26j in low energy regions (or near quark 
thresholds), where nonperturbative effects come into play. One thus defines 



D(-s) 



E 

0=0,2,4, 



y: c D {8,ii)(cb(ji)), 



(16) 



dimO=D 



where the Wilson coefficients C(s, /1) include short- distance effects and the operator 
(Ot>([i)) collect the long-distance, nonperturbative dynamical information at the arbi- 
trary separation scale fi. The operator of Dimension D = (D = 2) is the perturbative 
prediction (with mass correction). The D = 4 operator are linked to the gluon and quark 
condensates. Effective approaches together with the vacuum saturation assumption are 
used to compact the large number of dynamical D = 6 (D = 8) operator into one phe- 
nomenological operator (Oq) ((Os))- The nonperturbative addition to the Adler function 
reads then [|6[ 



^ of 27r 2 / 11 a q (s 





^GG 



18 71 I s< 

(m f q f q f ) 32vr 2 a s (s) 



27 7T 



E 



(m k q k q k ) 



(17) 



with the gluon condensate, ((a s /7i)GG), and the quark condensates, (rrifqfqf). The latter 
obey approximately the PCAC relations 



{m u + m d ) (uu + dd) 
m s (ss) 



n2 2 



where f n = (92.4 ± 0.26) MeV |2^] is the pion decay constant. The complete dimension 
D = 6 and D = 8 operator are parametrized phenomenologically using the vacuum ex- 
pectation values (Oq) and (0 8 ), respectively. Note that in zeroth order a s , i.e., neglecting 
running quark masses, nonperturbative dimensions do not contribute to the integral in 
Eq. (|TT|) . Thus in the formula presented in Eq. (|T8|) only the gluon and quark condensates 
contribute to R via the logarithmic s-dependence of the terms in first order a s . 

The total Adler D-function then reads as the sum of perturbative, mass and nonper- 
turbative contributions: 



D(s) = Dp(s) + D B 



D NP (s) 



(19) 



Uncertainties of the Theoretical Prediction 

Looking at Eq. (|T^) it is instructive to subdivide the discussion of theoretical uncertainties 
into three classes: 
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'i) The perturbative prediction. The estimation of theoretical errors of the perturbative 
series is strongly linked to its truncation at finite order in a s . Due to the incomplete 
resummation of higher order terms, a non- vanishing dependence on the choice of the 
renormalization scheme (RS) and the renormalization scale is left. Furthermore, one 
has to worry whether the missing four-loop order contribution d^(a s / 7r) 4 gives rise to 
large corrections to the series (0). On the other hand, these are problems to which 
any measurement of the strong coupling constant is confronted with, while their 
impact decreases with increasing energy scale. The error on the input parameter a s 
itself therefore reflects to some extent the theoretical uncertainty of the perturbative 
expansion in powers of a s . 

Let us use the following, intrinsically different a s determinations to benchmark 
our choice of its value and uncertainty. A very robust a s measurement is obtained 
from the global electroweak fit performed at the Z-boson mass where uncertainties 
from perturbative QCD are rather small. The value found is ct s (M|) = 0.120 ± 
0.003 0. A second precise a s measurement is obtained from the fit of the OPE 
to the hadronic width of the r and to spectral moments 0. The measurement is 
dominated by theoretical uncertainties, from perturbative origin. In order to ensure 
the reliability of the result, i.e., the applicability of QCD at the r mass scale, spectral 
moments were fitted analogously to the analysis presented in the following section. 
The nonperturbative contribution was found to be lower than 1%. Additional tests 
in which the mass scale was reduced down to 1 GeV proved the excellent stability 
of the a s determination. The value reported by the ALEPH Collaboration |J is 
a s (M|) = 0.1202 ± 0.0026. A third, again different approach is employed when 



using lattice calculations fixed at bb states to adjust a s . The value given in Ref. [pq| 
is a s (M|) = 0.117 ±0.003. 

The consistency of the above values using quite different approaches at various 
mass scales is remarkable and supports QCD as the theory of strong interactions. 
To be conservative, we choose a s (M|) = 0.1200 ± 0.0045 as central value for the 
evaluation of the perturbative contribution to the Adler D-function. 

Even if it is in principle contained in the uncertainty of Aa s (Mf) = 0.0045, we 
furthermore add the total difference between the results obtained using FOPT CI 
and those from FOPT as systematic error. 

(ii) The quark mass correction. Since a theoretical evaluation of the integral (|5]) is only 
applied far from quark thresholds, quark mass corrections -D mass are small. Without 
loss of precision, we take half of the total correction as systematic uncertainty, i.e., 
we add D mass ± D mass /2 in Eq. QT^) . 

(Hi) The nonperturbative contribution. In order to detach the measurement from theoret- 
ical constraints on the nonperturbative parameters of the OPE, we fit the dominant 
dimension D = 4, 6, 8 terms by means of weighted integrals over the total e + e~ low 
energy cross section. Again, without loss of precision, we take the whole nonpertur- 
bative correction as systematic uncertainty, i.e., we add D^p ± -Dnp in Eq. ([T9|). 

Another sources of tiny uncertainties included are the errors on the Z-boson and the top 
quark masses. 
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5 Spectral Moments 



Constraints on the nonperturbative contributions to R(s) from theory alone are scarce. 
It is therefore advisable to benefit from the information provided by the explicit shape 
of the hadronic width as a function of s in order to determine the magnitude of the OPE 
power terms at low energy. We consequently define the following spectral moments 



, 7 ds /_ s\Vs v ' 



R kl (s )= / - 1-- - R(s) , (20) 
J s \ s J \s / 

where the factor (1 — s/so) k squeezes the integrand at the crossing of the positive real axis 
where the validity of the OPE is questioned. Its counterpart (s/so) 1 projects on higher 
energies. The new spectral information is used to fit simultaneously the phenomenological 
operators ((a s /ir)GG), (Oq) and (0$), a procedure which requires at least 4 — better 5 - 
input variables considering the intrinsic strong correlations between the moments which 
are reinforced by the experimental correlations and by the correlations from theoretical 
uncertainties. 

To predict theoretically the moments, one uses the virtue of Cauchy's theorem and 
the analyticity of the correlator II(s), since a direct evaluation of the integral (^) in the 
framework of perturbation theory (and even OPE) is not possible. With the relation (|9]), 
Eq. ( ^0|) becomes 

R kl (s ) = Gni i-(l- -t (-) 11(a) , (21) 

\s\=so 

and with the definition ( |i~0D of the Adler D-function one further obtains after integration 
by parts 



1 r ds 

: f 

\s\=SQ 

where D(s) is obtained from Eq. (|19|) . 



sl=sn L 



D(s) , (22) 



6 Data analysis and Determination of the Moments 

Due to the suppression of nonperturbative contributions in powers of the energy scale s, 
the critical domain where nonperturbative effects may give residual contributions to R(s) 
is the low-energy regime with three active flavours. We thus choose the energy scale So 
of the fit equal to the energy scale, where the theoretical evaluation of R(s) shall start. 
As demonstrated in isovector vector r decays 0, the scale = M T is an appropriate 
scale where nonperturbative effects are present, but essentially controlled by the OPE. 
In our case we manipulate isovector and isoscalar vector hadronic final states, i.e., more 
inclusive data, and might expect smaller nonperturbative contributions. We therefore set 
the energy cut to ^/so = 1.8 GeV. Up to this energy, R(s) is obtained from the sum of 
the hadronic cross sections exclusively measured in the occurring final states. 
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Moments (k,l) 


Data 


A(stat. + sys.) 


A(theo.) 


(2,0) 


0.770 


0.013 


0.024 


(2,1) 


0.204 


0.005 


0.012 


(2,2) 


0.073 


0.003 


0.001 


(2,3) 


0.035 


0.002 


0.001 


(2,4) 


0.020 


0.001 


0.001 




Moments (k, I) 


(1,0) 


(1,1) (1,2) 


(1,3) (1,4) 


(1,0) 


1 


0.98 0.75 


0.60 0.49 


(1,1) 




1 0.79 


0.67 0.57 


(1,2) 




1 


0.97 0.91 


(1,3) 






1 0.99 


(1,4) 






1 



Table 1: Spectral moments measured at ^ l /so = 1.8 GeV with experimental and theoreti- 
cal errors. Below, the corresponding correlation matrix containing the quadratic sum of 
experimental and theoretical covariances. 

The data analysis follows exactly the line of Ref. l n addition to the e + e~ annihi- 
lation data we use spectral functions from r decays into two- and four final state pions 
measured by the ALEPH Collaboration p9[ . Extensive studies have been performed in 
Ref. in order to bound unmeasured modes, such as some KK or the 7r + 7r~47r° final 
states, via isospin constraints. We bring attention to the straightforward and statisti- 
cally well-defined averaging procedure and error propagation used in this paper as in the 
preceeding one, which takes into account full systematic correlations between the cross 
section measurements. All technical details concerning the data analysis and the integra- 
tion method used are found in Ref. (T| . 

The experimental determination of the spectral moments (p0"D is performed as the sum 
over the respective moments of all exclusively measured e + e~ final states (completed by r 
data). We chose the moments k — 2, I — 0, • • • , 4, in order to collect sufficient information 
to fit the three nonperturbative degrees of freedom. Neglecting the s-dependence of the 
Wilson coefficients in Eq. (|16"D, the respective nonperturbative power terms contribute to 
the following moments: the dimension D = 4 term contributes to / = 0, 1, the D = 6 term 
to / = 0, 1, 2 and the D = 8 term contributes to the I = 0, 1, 2, 3 moments. The I = 4 
moment receives no direct contribution from any of the considered power terms. However 
its use is not obsolete, since it constrains the power terms through its correlations to the 
other moments. Table [I] shows the measured moments together with their (statistical 
and systematic) experimental and theoretical errors. Additionally given is the sum of the 
experimental and theoretical correlation matrix as it is used in the fit. 

Using as input parameters a s (M|) = 0.1200 ± 0.0045, yielding for three flavours 
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Agl = (372 ±76) MeV), and for the mass of the strange quarkflm s (l GeV) = 0.230 GeV, 
while setting m u = m d = 0, the adjusted nonperturbative parameters are 

(0.037 ±0.019) GeV 4 , 
-(0.002 ± 0.003) GeV 6 , 

(0.002 ± 0.003) GeV 8 , (23) 

with % 2 /d.o.f = 1.5/2. The correlation coefficients between the fitted parameters are 
p(((a 9 MGG),(0 6 )) = -0.41, p(((a s /n)GG),(0 8 )) = 0.55 and p((0 6 >, (O s )) = -0.98. 
As a test of stability we have additionally fitted the k — 2, 1 = 0, . . . , 4 spectral moments 
at ^/so = 2.1. The difference between the results of this fit and Eq. fl23|) is included 
into the parameter errors given in Eq. (|23|). With these values, the corrections to R(sq) 
at v^o = 1-8 GeV amount to 0.16% from the strange quark mass and < 0.1% from the 
nonperturbative power terms. The gluon condensate can be compared to the standard 
value obtained from charmonium sum rules, ((a s /n)GG) = (0.017 ± 0.004) GeV 4 P5| , 
which lies below our value. However, another estimation using finite energy sum 
rule techniques on e + e~ data gives the value of ((a s /ir)GG) = (0.044 to^i) GeV 4 in 
agreement with the result (0). Fitting the moments when fixing the dimension D = 
6, 8 contributions reduces the gluon condensate to 0.010 ± 0.002. One may additionally 
compare the fitted dimension D = 6, 8 operator to the results obtained from the r vector 
spectral functions, keeping in mind that only the isovector amplitude contributes in this 
case and thus the more inclusive isoscalar plus isovector moments from e + e~ annihilation 
are expected to receive smaller nonperturbative contributions^. Reexpressing the results 
of Ref. |J in terms of the definition adopted in Eq. JT7|) we obtain (0 6 ) 7=1 = — (0.0042 ± 
0.0006) GeV 6 and (0 8 )i=i = (0.0062 ± 0.0007) GeV 8 . 

As a cross check of the spectral moment analysis, we fit the three nonperturbative 
parameters and a s . The theoretical error applied reduces essentially to the theoretical 
uncertainties of the QCD perturbative series estimated in Ref. to be Aa s (M|) = 0.0023 
at the t mass scale which is of the same magnitude as the scale ^/so = 1.8 GeV used here. 
The fit of the k = 2, 1 = 0, . . . , 4 moments yields a s (M|) = 0.1205 ± 0.0053 (x 2 /d.o.f = 
1.4/1), in agreement with the values from other analyses cited above. The values for 
the gluon condensate and the higher dimension contributions are consistent with those 
of Eq (^3|). Fitting a s (s ) and the gluon condensate at y/so = 2.1 GeV when fixing the 
dimension D = 6, 8 contributions at the values of Eq. (p3|) yields a s (M|) = 0.121 ± 0.011 

The mass of the strange quark used in this analysis takes one's bearings from the recent experimental 
determination using the hadronic width of r decays into strange final states, R T .s, reported by the 
ALEPH collaboration |o| to be m 3 (l GeV) = 235±H MeV. The error on this mass has no influence on 
the present analysis since, conservatively, 50% of the total mass contribution given in Eq. ( |l5| ) is taken 
as corresponding systematic uncertainty. 

5 The vacuum saturation hypothesis offers a relation between the dimension D — 6 contribution and 
the light quark condensates ||. Using the formulae ( p7| ) and (|l^) one has 

004 

(O e ) ~ -—naMim) 2 « -lO" 3 GeV . 



a* 

7T 



GG 

(Oe) 
(Os) 
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and ((a s /ir)GG) = (0.042 ± 0.003) GeV. The consistency of these results with those 
obtained at t/so = 1.8 GeV supports the stability of the OPE approach within the energy 
regime where it is applied. Varying the cm. energy so and the weights k, I of the moments 
used to fit the nonperturbative contributions gives a measure of the compatibility between 
the data and the OPE approach. Differences found between the fitted parameters are 
included as systematic uncertainties in the errors of the values fl2"3"|). 



7 Evaluation of Aa had (M|) and a™ d 

In order to get the most reliable central value of the perturbative R(s) prediction which 
enters the integrals ([5]) and ([?]), we use the whole set of formulae given in Ref. pi] , includ- 



ing mass corrections up to order a 2 s . Despite this theoretical precision, the uncertainties 
keep conservatively estimated as described in Section |3|. 

Since deeply nonperturbative phenomena are not predictable within the OPE ap- 
proach, we use experimental data to cover energy regions near quark thresholds. The 
low energy results for Aahad(Af|) and a^ ad from Ref. |]] including r data are taken for 
\fs < 1.8 GeV. The narrow u, <f), J/ip and T resonances are parametrized using relativis- 
tic Breit-Wigner formulae as described in Refs. [32 , [I]]. In addition, R measurements of 



the continuum contributions in the environment of the cc threshold are taken from the 
experiments specified in Ref. [[![]. The technical aspects of the integration over data points 
are also discussed in Ref. 0. No continuum data are available at bb threshold energies in 



the range of 10.6 GeV < E < 12 GeV. A recent analysis of the T system |3B| however 
showed that essentially within the perturbative approach of Ref. |2l| it is possible to pre- 
dict weighted integrals over the bb states. The uncertainty of this approach corresponds 
to an estimated error of Aa s (M|) = 0.008. We therefore use this uncertainty of a s for 
the QCD prediction at bb threshold energies. Included is a small uncertainty originating 
from scale ambiguities when matching the effective theories of four and five flavours. 

For the theoretical evaluation of the integrals (|5|), (0) (via Eq. (|TT|)), we use the fol- 
lowing variable settings: 

M z = (91.1867 ±0.0020) GeV @, 
a(0) = 1/137.036 , 
a s (M|) = 0.1200 ± 0.0045 , 

m u = m d = , 

m s (l GeV) = 0.230 GeV , 

^c{ m c) — 1-3 GeV , 

Tnb{mb) = 4.1 GeV , 

m t (m t ) = (175.6 ± 5.5) GeV [0 , 



and the values (1231) for the nonperturbative contributions. 

There are no errors assigned to the light quark masses, since the half of the total quark 
mass correction is taken as systematic uncertainty. The error on mt is needed in order to 
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-Emm - -Emax (GeV) 


Aa had (M|) x 10 4 
Data Theory 


a 


fl had x 1Q io 


Ami -1-8 


56.9±l.lW 


- 




636.49 ± 7.41 (D) 


1.8 - 3.700 


32.4 ±3.1 


24.50 ±0.33W 


2.5 


33.84 ±0.53 (T) 


V>(3770) 
3.700 - 5.000 


0.29±0.08 w 
15.8 ± 1.7 W 


17.06 ±0.58 


0.5 


0.17 ±0.02 (D) 
6.93 ±0.62 (D) 


5.000 - 10.500 




a-\ _i_ n /in(*) 


1.0 


7.43 ± 0.09 (T) 


T(4S, 10860, 11020) 
10.500 - 12.000 


0.38 ±0.10 
7.6 ±0.5 


8.19±0.32W 


0.4 


0.55 ±0.03 (T) 


12.000 - 40.000 


75.2 ±2.7 


77.96 ± 0.30W 


1.0 


1.64 ±0.02 (T) 


40.000 - oo 




41.98 ±0.22W 




0.16 ±0.00 (T) 


J/V>(1S,2S) 


9.68 ± 0.68W 






7.80 ± 0.46 (D) 


T(1S,2S,3S) 


0.98 ±0.15^ 






0.09 ±0.01 (D) 


Ami ~~ 00 


277.8 ± 2.2 cxp ± 1.4 thco 


695.1 ±7.5 cxp ±0.7 theo 



(*' Value used for final result of Aahad(Af|) (last line). 



Table 2: Contributions to A«had(^z) an< ^ ^° a |T d f rom the different energy regions. 
The "a" column gives the standard deviation between the experimental and theoretical 
evaluations of Aahad(^l)- "(D)/(T)" in the last column stands for evaluation using 
Data/ Theory. 

estimate the systematic uncertainty of the matching scale when turning from five to six 
flavours. Table || shows the experimental and theoretical evaluations of Aahad(^l) an d 
fl had f or reS p ec tive energy regimes. The upper star denotes the values used for the 
final summation of A«had(^z) given in the last line. At E ~ 3.8 GeV a non-resonant 
DD production might contribute to the continuum^. We therefore use experimental data 
to cover energies from 3.7-5.0 GeV. 

A 20% correlation is assumed between the analytic evaluations of the narrow reso- 
nances where in each case a Breit-Wigner formula is applied. The theoretical errors are 
by far dominated by uncertainties from a s and the difference FOPT/FOPTci- For in- 
stance, the first energy interval where theory is applied (E £ {1.8 — 3.7 GeV}) receives 



6 Such a contribution must be tiny since it is suppressed by its form factor and the (1 — 4A/£,/s) 3 / 2 
threshold behaviour of a pair of spin particles. 
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the error contributions 



A«had(M|) 



x 10 4 = 24.502 ± 0.332 



0.243 - Aa s 

0.223 — perturbative prediction 

0.036 — quark mass correction 

0.003 — nonperturbative parts 

< 0.001 - AM Z 



where the small contributions from quark masses and nonperturbative dimensions show 
that the perturbative QCD calculation is very solid here. Remember that only Ins- 
dependent nonperturbative terms contribute to R(s). Theoretical errors of different en- 
ergy regions are added linearily except the uncertainties at bb threshold that are (partly) 
from individual origin so that a 50% correlation to other energy regimes is estimated here. 
Looking at Table |2| one notices the remarkable agreement between experimental data and 
theoretical predictions of Aah a d(^f|) even in the cc quark threshold regions where strong 
oscillations occur. The experimental results of R{s) and the theoretical prediction are 
shown in Fig. |l|. The shaded bands depict the regions where data are used instead of 
theory to evaluate the respective integrals. Good agreement between data and QCD is 
found above 8 GeV, while at lower energies systematic deviations are observed. The R 



measurements in this region are essentially provided by the 772 |34[] and MARK I [35|| 
collaborations. MARK I data above 5 GeV lie systematically above the measurements of 
the Crystal Ball 



and MD1 [R7[ Collaborations as well as the QCD prediction. 



8 Results 



According to Table the combination of the theoretical and experimental evaluations of 
the integrals (H) and (^) yield the final results 



A« had (M|) 


= (277.8 ± 2.2 cxp ± 1.4 theo ) x 10" 4 


a~\Ml) 


= 128.923 ± 0.030 exp ± 0.019 theo 


had 
a fj, 


= (695.1 ± 7.5 cxp ± 0.7 theo ) x 10- 10 


af 1 


= (11 659 164.6 ± 7.5 cxp ± 4.1 thco ) x lO" 10 



(24) 



The total a^ M value includes an additional contribution from non-leading order hadronic 



vacuum polarization summarized in Refs. |38|, |T| to be a^ ad [(« / 7r) 3 ] 



-16.2 ± 4.0) x 



10 



~ 10 , where the error originates essentially from the uncertainty on the theoretical eval- 

| pE0[] . Fig. shows a com- 



uation of the light-by-light scattering type of diagrams 
pilation of published results for the hadronic contribution A«had(^z)- Some authors 
give the hadronic contribution for the five light quarks only and add the top quark 
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Figure 1: Inclusive hadronic cross section ratio in e + e~ annihilation versus the cm. 
energy y/s. Additionally shown is the QCD prediction of the continuum contribution as 
explained in the text. The shaded areas depict regions were experimental data are used 
for the evaluation of Aahad(^f) an d a )T d ^ n addition to the measured narrow resonance 
parameters. The exclusive e + e~ cross section measurements at low cm. energies are 
taken from DM1 , DM2, M2N, MSN, OL YA , CMD, ND and r data from ALEPH (see Ref 
for detailed information) . ^ 



t 1 1 r 
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O 
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~| 1 1 1 1 1 1 1 r 

Lynn et al. 87 
Eidelman, Jegerlehner 95 
Burkhardt, Pietrzyk 95 
Martin, Zeppenfeld 95 
Swartz 96 
Alemany et al. 97 
This work 



_i i i i i i i i_ 



_L 



_L 



_l I L_ 



270 



280 



290 



300 



310 



320 



Figure 2: Comparison of estimates of and Aah a d(^I/ 
Refs. P, 0,11, 11,11 1- 



Ac S M z) ( Xl0_4 ) 



The values are taken from 



contribution of Aa to p(^f) = — 0.6 x 1CT 4 separately. This has been corrected for in 
the figure. The present inclusion of theoretical evaluations yields an improvement in 
the Aahad(A^I) precision of a factor of two compared to our previous determination, 
a _1 (M|) = 128.878 ± 0.085 The improvement of the accuracy on a^ ad with respect to 
Ref. |l[ where we found a^ ad = (701.1 ± 9.3) x 10~ 10 is also significant. Other estimates 
of a!! ad can be found in Refs. [|32l Eq, fl6|, |47 



The analysis using the spectral moments and also the measurement of a s at the r mass 
scale showed concordantly that the OPE can be safely applied in order to predict inte- 
grals over inclusive spectra at relatively low energy scales. Nonperturbative contributions 
are indeed tiny and well below the envisaged precision. This approach leads to a large 
improvement in precision, as compared to the method used in Ref. |^3| where perturbative 



QCD was assumed to be valid only above 3 GeV. Also the value used for a s (M|) in the 
latter analysis was less precise than the current value. Finally, better and more complete 
experimental information, in particular at low energy, is used in the present work. Exper- 
imental data on R are still necessary at very low energies and at the cc quark production 
threshold. Uncontrolled nonperturbative effects spoil the hadronic spectra and make the 
OPE approach unreliable. Thus, the final results on both, Acthad(^z) an d a )T d ' are s ^ n ^ 
dominated by experimental uncertainties, where in particular the unmeasured 7r + 7r~47r° 
and KK7T7T low energy final states, which must be conservatively bound using isospin 
symmetry 0, as well as the two-pion final state in the latter case, contribute with large 
errors. Also more precise data of the cc continuum at threshold energies are needed. 

We use the improved precision on a(M|) to repeat the global electroweak fit in or- 
der to adjust the mass of the Standard Model Higgs boson, M H igg S . As electroweak and 
heavy flavour input parameters we use values measured by the LEP, SLD, CDF, DO, 
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CDHS, CHARM and CFFR collaborations, which have been collected and averaged in 
Ref . . The prediction of the Standard Model is obtained from the ZFITTER electroweak 
library ||48|| . The fit adjusts M%, M top and a(Mz) which are allowed to vary within their er- 



rors. Freely varying parameters are a s and M H ig gs . We obtain a s (M|) = 0.1198 ± 0.0031 



in agreement with the experimental value of 0.122 ± 0.006 |49| from the analyses of QCD 
observables in hadronic Z decays at LEP. The fitted Higgs boson mass is 129lg2 3 GeV 
with x 2 — 16.4/15, compared to 105lg2 2 GeV when using the previous value of a(M|) 
from Ref. An additional error of 50 GeV should be added to account for theoretical 
uncertainties |48|]. Doing so, we obtain an upper limit for M H igg S of 398 GeV at 95% CL. 

Fig. |3| depicts the variation of A% 2 as a function of the Higgs boson mass for the new 
and previously used values of a(M|) JE]. 

Table H shows the dominant uncertainties of the input values of the Standard Model fit 
expressed in terms of Asin 2 ^^. This work reduces the uncertainty of a(M|) on MHi ggs 
below the uncertainties from the experimental value of sin 2 ^^* or from the Standard 
Model, i.e., theoretical origin. 



9 Conclusions 

We have reevaluated the hadronic vacuum polarization contribution to the running of the 
QED fine structure constant, a(s), at s = M| and to the anomalous magnetic moment of 
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Experiment 
«(MJ) 



0.00023 



0.00023 



this work 



0.00010 



Theory 

^Higgs 



0.00018 
0.00014 



0.00160 [M^ 



60-1000 GeV] 



Fit 



129 



■38 s 



GeV 



Table 3: Dominant uncertainties of input values of the Standard Model electroweak fit 



the muon, a M . We employed perturbative and nonperturbative QCD in the framework of 
the Operator Product Expansion in order to extend the energy regime where theoretical 
predictions are reliable. In addition to the theory, we used data from e + e~ annihilation and 
r decays to cover low energies and quark thresholds. The extended theoretical approach 
reduces the uncertainty on Aahad(A^I) by more than a factor of two. Our results are 
Aa h ad(Mf ) = (277.8 ± 2.6) x 10~ 4 , propagating a _1 (0) to a~ l {Ml) = (128.923 ± 0.036), 
and a^ ad = (695.1 ± 7.5) x 10 10 which yields the Standard Model prediction a^ M = 
(11659164.6 ± 8.5) x 10~ 10 . The new value for a(M|) improves the constraint on the 
mass of the Standard Model Higgs boson to MHi ggs = 129lg2 3 GeV. 
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